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ABSTRACT 

We present the results of quantum calculations based on entirely ab initio 
methods for a variety of molecular processes and chemical reactions involving the 
LiHe + ionic polar molecule. With the aid of these calculations we derive accurate 
reaction rates and fitting expressions valid over a range of gas temperatures 
representative of the typical conditions of the pregalactic gas. With the help of a 
full chemical network, we then compute the evolution of the abundance of LiHe + 
as function of redshift in the early Universe. Finally, we compare the relative 
abundance of LiHe + with that of other polar cations formed in the same redshift 
interval. 

Subject headings: Early Universe - Astrochemistry - Molecular processes - ISM: 
molecules - Methods: numerical - evolution 
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Introduction 



The suggestion that chemical processes involving lithium atoms, as part of simple 
diatomic species formed with partner atoms and ions of expected significant abundances 
(i.e. H, He, H + and He + ), could play a role in the evolution of the early Universe has been 



put 



'orward several years ago (e.g. see 



Lepp fc Shull 



1983 



Stancil et al. 



1996 



Galli fc Palla 



19981 ). The formation of Li-bearing molecules results from the production of Li atoms a few 



minutes after the Big Bang and the concur rent formation of other light atoms like H, D 



and He (I Wagoner et al. 



1967 



Peebles 



19931 ). The fractional abundances of these elements 



are sensitive to the values of the baryon density and temperature of t he early Uniyerse 



and therefore imply specific physical constraints to their actual values ( jCyburt et al. 



2008|). 



In fact, as the Universe expanded its radiation temperature decreased and the atomic 
ions which had originated from the above elements gave rise to neutral atoms from ionic 
recombination with ambient electrons, thereby producing molecular species chiefly by 
radiative association. Neutral molecules include H2, HD, and LiH. 

Because of its large dipole moment and low ionization potential, LiH has been 
considered a potential candidate for inducing spatial and/or spec tral distort i ons of the 



Cosmic Background Radi ation (CBR) , as or igina 



observationally tested by 



Maoli et al- 



ly suggested 



^y 



Persson et al. 



Dubrovichl ( 119931 ) and 



(12010T ). It is, in fact, its 



(119941 ) and 

possible role as a molecular coolant of the primordial gas because of the efficient radiative 
decay down its manifold of rotovibration al levels that made its lik ely presence a very 



important issue in several earlier studies ( jBougleaux &: Galli 



19971 ). By the same token, 



other molecular species formed as cations after atomic ionization during gravitational 
collapse, e.g. LiH + , HeH + , HD + and LiHe + , could also play a similar role as additional 
molecular coolants. Furthermore their non-equilibrium level populations may also have left 
a possible signatures in protogalactic clouds, imprinting spatial or spectral distortions in 
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the CBR spectrum (ISchleicher et al.ll2008l ). 



In particular, in this work we analyse LiHe + via a series of quantal calculations to 
determine the relative role of several pathways to its format ion and destruction within the 



2011a. 



20121 ) we 



chemical network acting at low redshift. In previous paper (jBovino et al.l 
examined the photonic paths to the formation and destruction of LiHe + . In the present 
paper we extend that work here by adding the chemical paths presiding over its evolution, 
i.e. the chemical reaction with hydrogen: 



LiHe H 



H -> LitT 



He 



leading to its destruction, and the formation reaction from the lithium hydride cation: 



LiH+ + He LiHe + + H, 



(2) 



the phy sical characteristics of which are discussed extensively in our recent work on this 



system (ITacconi et al. 



20JJ). 



An additional path that we shall further consider in the prese nt study is that dri ven 



by the presence of a residual electron fraction in the early Universe (jStancil et al. 



1996| ), i.e. 



the fragmentation of the polar cation by the dissociative recombination (DR) process: 



LiHe + + e" -> LiHe* -)■ Li + He + hu 



a reaction that, from the present calculations, shall turn out to be fairly efficient. 



(3) 



In a previous study of the lithium chemistry in the early Universe (jBovino et al 



2011bl ). we found that the abundances of LiH and LiH + were of the order of about 10~ 17 
(number fraction), i.e. too diluted for direct detection with presently available tools. Hence 
it is of interest to see what the abundance would be for the LiHe + , a molecular partner not 



considered in previous studies because of lack of reliable data. To this goal, we shall avail 
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ourselves of the new, and accurate, cross sections and reaction rates obtained here with 
fully quantal methods. 

In the next section we will briefly summarize the computational findings for reactions 
(1), (2), and (3) and report the corresponding computed reaction rates at the relevant 
temperatures, while in section 3 we will discuss the relative role of the formation/destruction 
channel rates of LiHe + at low z, together with its fractional abundance compared to those 
of other molecules present in the early Universe. Finally, in section 4 we shall summarize 
our conclusions. 



2. The computed quantum rates 



As mentioned in the introduction, our previous study on the pho ton- induced evolution 



of LiHe + under early Universe conditions (IBovino et al. 



2011 



a, 



2012l ) . has shown efficient 



formation of the latter molecule by resonant radiative recombination but also competitive 
destruction by photodissociation at low redshifts. Hence, it was suggested by the ensuing 
modeling of its relative abundances that the LiHe + number fraction at z ~ 30-1 would be 
lower than 10 -18 by photon-induced paths only, i.e. much lower that the corresponding 
values obtained for LiH + (~ 10~ 14 ) and for HeH + (~ 10~ 9 ) under the same conditions. On 
the other hand, a more realistic picture should be obtained by the additional implementation 
to the evolution scheme of the chemical routes outlined in the previous section, in order to 
achieve a fuller comparison with the data which already exist for the other polar cations 
formed in the early Universe. 

The chemical reaction route for both processes: 



LiH+( 2 S + ) + He( 1 S) ^ LiHe+( 1 £+) + H( 2 S) 



(4) 



of formation (with a threshold) and of destruction (without threshold) have been studied in 
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detail earlier by us (ITacconi et al.ll20 



for the reactive paths (j Wernli et al. 



2|) us ing accurately computed potential energy surfaces 



20091 ). The quantum reactive cross sections were 



obtained using a Negative Imaginary Potent ial (NIP) approach w ithin a Coupled-States 



(CS) description of the quantum dy namics ( 



computati ons have been reported by 



Bovino et al. 



Tacconi et al 



201 Id ) . Full details of the 



reactions (jBovino et al. 



2010 



(120111 ) and in earlier papers on ionic 



201ld ): we therefore refer the interested reader to those 



papers for further details. In order to put the present results in a more specific context, 
we show in figure 1 the behaviour of the computed rates for the formation (left panel) 
and destruct ion (right panel) re actions given by eq. (4) and obtained from our recent 



calculations ( ITacconi et al. 



20121 ). 



It is interesting to note at the outset the contrasting beha viour between the fo rmation 



20121 ) and 



of the cation, which presents an energy barrier of ~ 0.05 eV (ITacconi et al. 
therefore becomes negligibly small at lower redshifts, and the same process when the formed 
LiH + is taken to be vibrationally excited (dashed line in the left panel). Since the amount 
of internal energy is sufficient to overcome the barrier to reaction, the formation process 
now shows a rather mild dependence on temperature and exhibits much larger formation 
rates around 10~ 10 cm 3 s -1 . 

Entirely different, however, is the behaviour of the exothermic destruction reaction, 
reported on the right panel of figure 1. In that case the reaction rates are largely 
independent of temperature but, even for LiHe + in its ground rotovibrational state, remain 
of the order of 3-4 xlO -10 cm 3 s _1 at the lower temperature values. It is therefore natural 
to expect that the different sizes of the two reaction channels will play an important role 
when generating the final abundances, as we shall discuss further in the next section. 

Another "chemical" process which we need to consider concerns the effects of the 
ambient electrons on the formed cationic molecule as described by reaction (3), i.e. the 
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destruction of LiHe + following the DR mechanism outlined there. The DR process has been 
studied before and correspondi ng data un der early Universe con ditions already exist for 



HeH + (IHaxton fc Greene 



2009|) and LiH+ flCurik fc Greene 



20071 ). No data, however, were 



available for the present system. 

Therefore we have analysed the electron-assisted destruction of LiHe + , with the 
ensuing pr oduction of Li an d He atoms (the details of the calculations will be presented 



elsewhere (jCur'ik et al 



201 21). Here, it suffic es to say that the recombinative process, similar 



to the case of LiH + (jCur'ik fc Greend 120071 ) is not a direct curve-crossing process with a 
dissociative, neutral potential but rather an indirect process which involves Rydberg states 
of the lithium atom perturbed by the presence of He. The final destruction is therefore 
caused by vibrational Feshbach resonances within the (Li + He)* collisional complex which 
lead to th e efficient break-u p between a radiatively stabilized lithium and a neutral helium 



atom (see 



Cur'ik et al 



2012 



for further details). 



The general behaviour of the computed rates as a function of the temperature of the 
gas T g (assumed equal to the electron temperature), is given by the calculations shown in 
Figure 2. The figure clearly indicates the marked dependence of such rates on the internal 
energy of the initial target, while the definition of a global equilibrium temperature is the 
one which can be used within the calculation of the relative abundances, as we shall report 
in the following section. 



3. The chemical network and the evolutionary modelling 



The evolution of the pregalactic gas is usually considered within the framework of 
a Friedmann cosmological model and the cosmological abundances of the main atomic 
components are taken from standard Big Bang nucleosynthesis results ( jSmith et al.lll993l ). 
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The numerical values of the cosmological 



from WMAP5 data (Komatsu et a 



Galli fc Pallal fll998f ): 



Coppola et al 



2009). 



parameters used in the calculation are obtained 
For additional deta ils on the model, see 



(2011) 



Bovino et al 



( 1201 lbh . 



In order to obtain the abundances of LiHe + which evolve from the network of the 
relevant processes, a set of differential coupled rate equations of the form: 

dm 



—jjT = aform^fc - a dest rii + ... 



(5) 



has been solved. In eq. (5), «f orm and o^est are the formation and destruction rates of 
the species under discussion, and is the number density of the reactant species i. The 
rate coeffici ents involved i n the lithium chemistry network are basically the same as those 



reported in 



Bovino et al. 



( 120121 ) . with the exception of the reactions involving LiHe + 
molecule which now come from the present quantum calculations. Table 1 shows the 
adopted reactions and the fitting formulae of the rate coefficients. The equations governing 
the temperature and redshift evolution s have been reported in earlier papers (see e.g 



Galli fc Palla 



1998 



Bovino et al.ll2011bl ) and will not be repeated here. 



4. LiHe + evolution and final abundances 

As mentioned in the previous sections, it is instructive to model, from the present 
data, plus the network of rate equations outlined before, the way in which the production 
of LiHe + molecules evolves, as a function of redshift, in its interplay with the coexisting 
destruction channels. 

The complex evolution of all the considered species is shown in Figure 3, where the 
labels pi indicate production channels, whereas the labels di identify the destruction channels 
(same notation as in Table 1). The production of LiHe + is dominated at all redshifts by 
radiative association of Li + with He, both spontaneous and stimulated (reactions p\ and 
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P2, respectively). The LiH + channel is never effective, because of a threshold at energies 
corresponding to ~ 40 K in the reaction of LiH + (v = 0) with He (reaction p 3 ). The reaction 
of LiH + (t> = 1) with He (reaction p^), although characterized by a much larger rate than 
its v = counterpart, is also not important because of the rapid decay of the v = 1 level 
population of LiH + at low redhifts. The destruction of LiHe + is due to photodissociation 
for z > 60 (reaction d\) and to collisions with H for z < 60 (reaction d 3 ), with a significant 
contribution from dissociative recombination (reaction d-z). 

The above data are considered reliable down to z values of about z = 10 since at 
that stage the first stars are formed and reionization occurs, thereby invalidating the 
assumptions at the basis of the chemical network. On the whole, however, we can say that 
the close competition between efficient formation paths, and nearly as efficient destruction 
reactions, indicates for the present system a limitation to its final abundances in the region 
of redshift of interest. 

The consequences of the interplay between the different conflicting rates discussed 
above could be more clearly seen in Figure 4 showing the different abundances of various 
molecular cations on a log/log scale as a function of redshift. Most of the data for species 
other than LiHe + have been presented in earlier work, while the evolution of LiHe + is 
obtained from the new, accurate quantum calculations of the present study. 

Due to the concurrent destruction paths, it is clear from that comparison that, 
at redshift values between about 20 and 10, the fractional abundance of the LiHe + 
ions never increases beyond ab out 10 -22 number fractions, well below those of HeH + 



and LiH + (IBovino et al. 



2011b 



|cj). While the former molecule clearly remains the most 
interesting candidate for experimental observation, LiH + is still a borderline case that 
provides a challenging option while LiHe + appears still below the sensitivity of current 
instrumentations. 



5. Conclusions 



In the work reported in this paper we have analysed in greater detail the molecular 
processes which involve a thus far poorly studied molecular cation, the LiHe + ( 1 S + ), 
and which deal with its possible formation and destruction in the pregalactic gas 
through an extensive network of photon-induced and chemically driven processes. In 
the range of redshifts of interest we have carried out quantum calculations of chemical 
formation/destruction reactions which have never been considered before from realistic 
computational models (see reactions (1) to (3) of section 2). 

Our results clearly show the close competition between production and destruction 
pathways which turn out to have very similar efficiencies. Hence, the estimated abundances 
for LiHe + , within z values from about 30 and 10, do not increase as highly as those found 
for LiH+ and HeH+. 
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Fig. I. — Computed quantum rates (from Tacconi et al. 2011a) of the LiHe + formation 
reaction {left panel) and of its destruction by hydrogen {right panel) down to the low redshifts 
of early Universe environment. 
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Fig. 2. — Computed dissociative recombination (DR) rates for LiHe + as a function of the 
gas (electron) temperature. Different vibrational levels of initial LiHe + are shown, while the 
solid black curve represents the case where a temperature balance exists between the internal 
energy of the ion and that of the electrons. 
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Fig. 3. — Computed production/destruction channel evolutions for LiHe + as a function of 
the redshift z. The present data for LiHe + include the corrections to the results by Bovino 
et al. (2011a) presented by Bovino et al. (2012). 
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Fig. 4. — Computed abundances (number fractions) of molecular cations in the early Uni- 
verse as a function of the redshift z, obtained from the rates reported in Figure 3. 
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Table 1: Reactions considered in the present study. Fitting expression as function of the 
gas (T g ) or radiation (T r ) temperatures are accurate to better than 10% in the range 1 K < 
T g , T r < 10 3 K. 





reaction 


fitted rate (cm 3 s 1 or s 1 ) 


ref. 


PI 


Li+ + He -> LiHe+ + hv (sp.) 


logfe = -21.87 - 0.204 logT 9 - 0.153(log T g ) 2 + 0.233(log T g ) 3 - 0.144(log T g ) 4 + 0.0148(log T g ) 5 


(a) 


VI 


Li+ + He -> LiHe+ + hv (st.) 


log(fc/T r ) = -22.23 - 0.600 log T g - 0.086(log T g ) 2 - 0.082(log T g ) 3 + 0.042(log T 3 ) 4 - 0.0073(log T g ) 5 


(a) 


P3 


LiH+(i> = 0) + He -t- LiHe+ + H 


fe = (5.047 x 10~ 12 - 2.723 X 10~ 14 T g + 7.775 X U)- 17 T^ - 1.954 x 10 _20 T|) cxp(-441/T g ) 


(b) 


P4 


LiH+(i> = 1) + He -+ LiHc+ + H 


logfe = -9.96 - 0.648 log T g - 2.743(log T g ) 2 + 6.092(log T g ) 3 - 9.297(log T g ) 4 + 8.397(log T g ) 5 








-4.420(logT 3 ) 6 + 1.328(logT g ) 7 - 0.210(log T g f + 0.013(log T g f 


(6) 


di 


LiHe+ + hv -> Li+ + He 


k = 21 T°- 25 cxp(-5460/T r ) 


(a) 


d 2 


LiHe+ + e -> Li + He 


logfe = -5.254 - 0.452 log T g - 0.791(log T g ) 2 + 0.463(log T g f - 0.076(log T g ) 4 


(c) 


d 3 


LiHc+ + H -> LiH + + He 


log fe = 5.780 - 2.274(log T g ) + 0.154(log T g ) 2 + 0.0639(log T g ) 3 


(a ) 



(a) Bovino, Tacconi, Gianturco (2012); (6) Tacconi, Bovino, Gianturco (2012); (c) Curfk ct al. (2012, in preparation). 



